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How mixed can one component of a bi-partite system be initially and still become entangled through inter- 
action with a thermalized partner? We address this question here. In particular, we consider the question of how 
mixed a two-level system and a field mode may be such that free entanglement arises in the course of the time 
evolution according to a laynes-Cummings type interaction. We investigate the situation for which the two-level 
system is initially in mixed state taken from a one-parameter set, whereas the field has been prepared in an ar- 
bitrary thermal state. Depending on the particular choice for the initial state and the initial temperature of the 
quantised field mode, three cases can be distinguished: (i) free entanglement will be created immediately, (ii) 
free entanglement will be generated, but only at a later time different from zero, (iii) the partial transpose of the 
joint state remains positive at all times. It will be demonstrated that increasing the initial temperature of the field 
mode may cause the joint state to become distillable during the time evolution, in contrast to a non-distillable 
state at lower initial temperatures. We further assess the generated entanglement quantitatively, by evaluating 
the logarithmic negativity numerically, and by providing an analytical upper bound. 



I. INTRODUCTION 

Interacting quantum systems typically develop correlations 
in the course of their dynamics, either of classical or of gen- 
uinely quantum nature [U|-H]. If two systems are initially in a 
pure product state which does not decouple from the interac- 
tion Hamiltonian, entanglement will always occur. Depending 
on the particular type of the interaction and the choice of the 
initial state, it may well be that the degree of entanglement 
being generated is very small indeed for a rather long time pe- 
riod |Q] . Entanglement is nevertheless an unavoidable byprod- 
uct of the time evolution. If the joint system is initially in a 
product state, but with one or both subsystems being mixed, 
then it is not necessarily clear that entanglement will arise at 
any point in the dynamics. For finite-dimensional systems it 
is obvious that whenever the joint system is suitably mixed 
initially, then the joint state will stay separable at all times. 
For infinite-dimensional systems, however, a more complex 
behaviour may be anticipated. Intuitively, one may still ex- 
pect that if both systems are initially very mixed, then the bi- 
partite state should not develop entanglement at any time. In 
the case that one system is in an arbitrarily mixed state and 
the other is pure, then it is conceivable that entanglement is 
always generated. This has in fact been demonstrated in the 
case of a Jaynes-Cummings-type interaction [|5|-|7p of a two- 
level system coupled to a field mode and in the Caldeira- 
Leggett model for quantum Brownian motion, where a single 
field mode is considered that is linearly coupled to many other 
field modes initially prepared in some thermal state [|^] . 

The Jaynes-Cummings model is also particularly suitable 
to assess the situation when initially both systems are initially 
in a mixed state. This model describes in the most basic way 
the interaction between light and matter: it models the dy- 
namics under linear coupling of a two-dimensional system to 
a near resonant quantised single field mode with a rotating- 
wave Hamiltonian. The Jaynes-Cummings model can be an- 
alytically solved in the sense that the time evolution can be 
explicitly evaluated for all times. Hence, one may investigate 



the trade-off between the degree of mixing of the two-level 
system and the mixing of the field mode. 
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FIG. 1. This figure depicts the maximal degree of entanglement 
that arises in the joint state of a two-level system coupled to a sin- 
gle field mode in the Jaynes-Cummings model. The field mode 
is initially in a thermal state, the two-level system is prepared in 
A|e)(e| + (1 — X)\g)(g\, A 6 [0, 1]. Zero temperature corresponds 
to a = 0, and the limit of infinite temperature is reflected by the 
limit a —> 1. The degree of entanglement is measured in terms of 
the logarithmic negativity. The shaded area represents the situation 
where the log-negativity is zero for all times, meaning that the joint 
state has a positive partial transpose at all times. In this figure, a is 
related to the initial inverse temperature (3 of the field mode accord- 
ing to a = m/(m + 1), where m = (exp(/3) — is the mean 
photon number. 

In our analysis, we consider the case where the field mode 
is initially in a thermal state. The set of initial states of the 
two-level system, from now on labeled as A, will be given by 

p A = A|e><e| + (1 - A)b><ff|, Ae[0,l], (1) 
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where |e) and \g) denote the state vectors corresponding to the 
excited and the ground state, respectively. We will analyse the 
entanglement properties of the joint system at all times, for all 
initial temperatures of the field mode and all A G [0, 1]. It will 
turn out that depending on the actual parameters in the model, 
three regimes can be distinguished: the regime where (free) 
entanglement is generated immediately, the regime where the 
joint state will be entangled at some later time t > 0, and the 
regime where the joint state has a positive partial transpose 
for all times. The regime where entanglement will immedi- 
ately be created is by no means restricted to the case where 
the two-level system is initially pure. At arbitrarily large tem- 
perature of the field mode one may still have a very mixed 
initial state of the two-level system - yet different from the 
maximally mixed state - in order to obtain an entangled state 
for all times. 



II. JAYNES-CUMMINGS MODEL FOR ATOM-FIELD 
INTERACTION 



The Jaynes Cummings interaction [gjjg] of the near reso- 
nant interaction with a single mode quantised field with anni- 
hilation and creation operators a and a*, respectively, can be 
written as 



H^t = 7(|e)(ff|o+|s)(e|at) 



(2) 



For simplicity, we set h = 1 and assume the frequency of the 
field mode to be unity as well, which does not imply a restric- 
tion of generality for our argument. The field mode, associ- 
ated with the label F, is initially prepared in a thermal state 
with respect to the inverse temperature > 0, i.e, 



exp(-/3H F ) ^ 



(3) 



n=0 



where Hp '■= a^a + 1/2 is the Hamiltonian of F and 



Pn 



(1 + m) n+1 ' 



(4) 



m is the mean photon number at the inverse temperature (3, 
and is given by 



m := (exp(/3) — 1) 1 



(5) 



The joint initial state ps®PF evolves in the interaction picture 
according to 



pit) = u(t) {p s ®PF)u{ty, t>o 



(6) 



where U(t) := exp(—iH int t). The formal time evolution 
can be made explicit, as the Jaynes-Cummings model can be 
analytically solved. The state pit") at a time t > can be 
evaluated using the standard techniques to solve the Jaynes 



Cummings model [R]. For brevity, we will write p instead of 
pit) from now on. With the above notation, one arrives at 



^PnPn, 



(7) 



where 



A 



c n \e,n){e,n\+s n \g,n + l){g,n + \\ (8) 
+ ic n s n (\g,n + l)(e,n\ - \e,n)(g,n+ 1|) 
+ (! - sl_x\e,n-l){e,n + c 2 n _ 1 \g,n)(g,n\ 

+ ic n -is n -i (|e,n- l){g,n\ - \g,n){e,n- 1|) , 

if n > 0, and 

Po=A[^|e > 0>(e,0| + ag|s,l>( ff> l| (9) 

- ic s \e, 0)(g, 1| + ic s \g, l)(e,0| 
+ (1-X)\g,0)(g,0\. 



The symbol 



Lo n :=2 7 (n+l) 1 / 2 



(10) 



has been used for the Rabi frequency, and s n and c„ are 
defined as c„ := cos(w„t/2), and s n := sin(w n i/2) for 
n = 0, 1, 2, .... We will set 7 = 1, again for reasons of clarity, 
which simply corresponds to a rescaling of the units of time. 



m. ENTANGLEMENT PROPERTIES AND PARTIAL 
TRANSPOSE 

In order to assess the entanglement properties of the joint 
state p at a later time t > we investigate the spectrum of the 
partial transpose p T of p. This partial transpose takes the form 
of a direct sum, such that the spectrum of p T can be evalu- 
ated from the spectrum of each block. In order to simplify the 
notation, we set 



p n = (1 - a)a n , 



(11) 



with a e [0, 1). The partial transpose is then explicitly given 
by 

OO 

p v = (1 - a) a"" 1 [(«^A + a 2 4(l - A))|e, n)(e, n\ 

n=0 

+ iiXa - (1 - A)a 2 ) 

x c n _ 1 s n _i(|e,n+ l){g,n\ - \g,n)(e,n + 1|) 



+ (4-i A + (! - X)c n -xOi)\g,n){g,n\ 



(12) 



From this equation, the spectrum er(p r ) of p T can be read off: 
it is given by 



(J (p r ) = {l -a>Ao,A , A+, A x ,...} , 



(13) 
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where A+,A n , n = 0,1,2,..., with A+ > A„ are the two 
eigenvalues of the 2 x 2-matrices 



A<n 

c: 



C n 
B„ 



A)), 



with entries 

A n := (l-a)a n+ \c 2 n+1 \ + as 2 n+1 (l 
B n := (1 - a)a"- 1 (o4_ 1 (l - A) + s^A), 
C„ := i(l — a)a n (\ — a(l — A))c„s„ 

if n > 0, and 



A 
Bo 

C 



= (1 - a)a(cf\ + as\(l - A)), 

= (l-a)(l-A), 

= i(l — a) (A — a(l — A))cqSo- 



(14) 



(15) 
(16) 
(17) 



(18) 
(19) 
(20) 



The remaining task is to identify those values for a and A for 
which there exists a t > such that at least one A~ is nega- 
tive. Due to the form of M n it is obvious that A+ > for all 
n = 0, 1, 2, .... Hence, it is sufficient to consider the positiv- 
ity of det[M„]. We first consider the case that n > 0. The 
numbers 



i 1 ' 2 , (n+l) 1 / 2 , (n-1) 1 



/2 



(21) 



are not in a rational relationship to each other. Therefore, 
each contribution to det[M n ] may be individually minimised. 
Hence, there always exists a time t > such that M n is given 
by Eq. Q with 



where 



A n = (l-a)a n+1 f(a,X), 
B n = (l-a)a n - 1 f{a,X), 
C n = »(1 - a)a n (\ - a(l - A))/2, 



/(a, A) := min(A, a(l — A)). 



(22) 
(23) 
(24) 



(25) 



This choice of t corresponds to the sitation where det[M ra ] 
takes its minimal value on the interval t E (0, oo). This means 
that there exists a time t > with det[M„] < if and only if 



/(«,A) 2 -(A-a(l-A)) 2 /4<0. 



(26) 



In the case n — 0, then one arrives with an analogous argu- 
ment at the statement that there exists a time t > such that 
det[M ] < Oif and only if 



af{a, A)(l - A) - (A - a(l - A)) 2 /4 < 0. 



(27) 



Eqs. (^) and ( j2~7j ) provide the desired relationship between a 
and A. Three cases may now be distinguished: 

(i) The values of A 6 [0, 1] and a € [0, 1) are such that 
neither Eq. ( 26 ) nor Eq. ( p7[ ) are satisfied. This means 
that the joint state of the two-level system and the field 
mode has a positive partial transpose for all times (i.e., 
it is a PPT state). 



(ii) The values of A and a are such that Eq. (p_7|) is satisfied, 
but Eq. (|6|) is not satisfied. This means that there exists 
a time t > such that the joint state p has a non-positive 
partial transpose (i.e., it is a NPPT state). This time can 
easily be evaluated: If 



A < a(l — A), 



(28) 



then it is given by the time t > that achieves cqSq = 
1/2 and ci = 1. If A > a(l - A), then c s = 1/2 and 
Si = 1 have to be simultaneously satisfied. 



(iii) The values of A and a are such that Eq. ( |2q ) is fulfilled. 
Then, again, the joint state will become a NPPT state. 
However, now the situation is different compared to 
case (ii): for any time t > one can find an n = 1,2,... 
such that det[A/ n ] < 0. This means that the joint state p 
becomes an NPPT state instantaneously, or, more pre- 
cisely, p is an NPPT state for all times t G (0, oo). 

It is important to note that whenever the state is a NPPT 
state, it is also a free entangled one. This can be seen when 
considering 



(1 <g> 7T n )p(l (g) 7T n ) 



fom = 1, 2, where 



n+l 



(29) 



(30) 



This is for each n a projection onto a 2 x 3-dimensional sub- 
space. By investigating the block diagonal structure of p r , 
one finds that the partial transpose of (1 ® n n )p(l ® ir n ) has 
again a direct sum structure. It follows from direct investi- 
gation of this projection that p T is non-positive if and only if 
there exists an n such that (1 (8 7r„)p(l ® 7r„) is non-positive. 
Hence, whenever one concludes that p is NPPT, then it is also 
free entangled, as for 2 x 3-systems being NPPT and being 
distillable are equivalent ||] . 

The situation that has just been described is depicted in 
Figs. 1 and 2. In Fig. 1, the white area corresponds to the 
situation where for some time t > a free entangled state is 
reached, which is case (ii). The shaded area is assigned case 
(i), the PPT case. The lines in Fig. 2 correspond to equality in 
Eq. (|26|): then free entanglement is created immediately, as in 
case (iii). 

It is interesting to note that even in the limit a — ► 1, which 
physically is the limit of infinitely large temperature, encoun- 
ters a PPT state. No matter how large the temperature, there 
always exists a value A £ (0, 1) such that the joint state be- 
comes entangled. The intersection of the parametrised curve 
of equality in Eq. ( p6| ) with a = 1 associated with 'infinite 
temperature' can be easily evaluated. One finds two intersec- 
tion points at A = 3/4 and A = 1/4. The intersection points 
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of the curve of equality in Eq. ( |27| ) with a — 1 are given by 
A = 1/2 — y/2/A = 0.14645 and A = 3/4. That is, whenever 
A is larger than 3/4 or smaller than 1/2 — \/2/4, then one 
can be sure that distillable entanglement will be generated. To 
this extent, both the two-level system and the field mode may 
therefore be mixed initially. If A = 1 /2, which corresponds 
to the maximally mixed state for the two-level system, there 
exists a temperature for which there is no distillable entangle- 
ment left. Another interesting situation arises when A = 1/10 
is chosen: Then, at small temperatures, distillable entangle- 
ment will be generated. As the temperature increases, that is, 
with increasing a, no free entanglement will occur. But then 
again, at very large temperatures free entanglement will again 
arise. 




FIG. 2. This figure depicts the boundary between the region 
where free entanglement is created immediately and the region for 
which there exists a t > such that the joint state is PPT at time t . 



From the above considerations, one can also find an up- 
per bound for the distillable entanglement with respect to lo- 
cal operations with classical communication that can maxi- 
mally be achieved in the course of the time evolution. The 
log-negativity is such an upper bound for the distillable en- 
tanglement. The negativity [p|Jlo|] of a state p is defined as 
||p r ||i — 1, where ||.||i denotes the trace-norm; it has been 
shown to be an entanglement monotone [foiO]. The log- 
negativity is nothing but N(p) := log 2 ||/? r ||i; it is an eas- 
ily computable upper bound for the distillable entanglement 
JlC|]. Fig. 1 depicts the maximal log-negativity with respect 
to all times from numerical investigations. Note that around 
A = 1/10, the counterintuitive situation occurs that the max- 
imal degree of entanglement that is achieved increases with 
the initial temperature of the field mode. 

An upper bound for the infimum N(p) :— log ||/? r ||i over 
all times can be evaluated from the above explicit expression 



for p T . First, note that 



= 1 



2£l A «l- 

n=0 



(31) 



An upper bound for the log-negativity and therefore for the 
distillable entanglement can then be evaluated from the ex- 
tremal sitations in each block of p T . One arrives at 



N(p) < log 2 (l + 2 min(0, \p~\) + 2min(0, |)) 
where pT and v~ are the smallest eigenvalues of 

af(a,X) (A-a(l-A))/2 



(1-a) 



(A-a(l-A))/2 



1-A 



and 



a 2 f(a,X) a{\-a{l- A))/2 \ 
a(A-a(l- A))/2 f(a, A) J 



(32) 



(33) 



(34) 



respectively. A particularly simple situation is the case that 
a = 1, corresponding to 'infinite temperature' of the field 
mode. Then one finds that 



f log 2 (2-4A), for AG [0,1/4], 
N(p) < { log 2 (4A-2), for A € [3/4,1], 
otherwise. 



(35) 



IV. CONCLUSIONS AND OUTLOOK 

In the main part of this paper we have investigated the prop- 
erties of the dynamically emerging entanglement in the one- 
photon Jaynes-Cummings model. We have found that three 
different regimes occur, depending on the actual initial joint 
product state: there may be free entanglement at all times, 
and only at a later time. Also, there is the case that the 
partial transpose of the joint state stays positive at all times, 
such that the state does not contain distillable entanglement 
That is, the state is then either separable or bound entan- 
gled. One should, however, keep in mind that for bi-partite 
systems, where at least one part is infinite-dimensional (such 
as the system considered in this paper), some care is needed 
when speaking about the separability of a state: the set of non- 
separable states is then trace-norm dense in the set of all states 
JT2] ] (see also Refs. [jT|]): this means that for any separa- 
ble state there exists an entangled state arbitrarily close to the 
original state as measured by the trace-norm, a state that leads 
consequently to arbitrarily similar probability distributions in 
quantum hypothesis testing distinguishing the two states. In 
contrast, there are clearly entangled neighbourhoods of entan- 
gled states. 

One - possibly not very surprising - principal observation 
is that the point corresponding to A = 1/2 and a = 1 in 
Fig. 1 is contained in the shaded region. This means that 
there always exists a temperature such that if both the two- 
level system S and the field mode F are initially in thermal 
states corresponding to this temperature, then the composite 
system will certainly not develop free entanglement. The spe- 
cific shape of the shaded area in Fig. 1 is clearly a property of 
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the considered one-photon Jaynes-Cummings model. By ap- 
propriately choosing the interaction Hamiltonian - for exam- 
ple by introducing a fc-photon Jaynes-Cummings interaction 
- one may obviously be in the position to change the shape 
of this area, and even to make its area very small. If one al- 
lows for a time dependence of the interaction, one can also 
find Hamiltonians such that for all initial temperatures of both 
parts of a bi-partite system, there will be a time t > such that 
the joint system will eventually become free entangled. This 
leads to the general question whether this connection is un- 
avoidable: there may even be bi-partite infinite-dimensional 
systems with the property that (i) the thermal states exist for 
all temperatures for both parts of the system and (ii) with the 
property that the composite system becomes entangled for all 
initial temperatures for a fixed time-independent interaction. 
The existence of such systems may not appear particularly 
plausible. A complete answer to this question may however 
lead to a more thorough understanding of the exact link be- 
tween the initial mixedness of the state of a composite quan- 
tum system and the capability of producing entanglement in 
the course of the time evolution. 



V. ACKNOWLEDGEMENTS 



This work has been supported by the European Union 
Networks EQUIP, QUEST, and QUBITS, the EPSRC, the 
DFG (SchwerpunktprogrammQIV), and the A.-v.-Humboldt- 
Foundation (Feodor-Lynen-Grants of JE and SS). 



[1] K. Zyczkowski, P. Horodecki, M. Horodecki, and R. 

Horodecki, P hys. Rev. A 65, 012 101 (2002); S. Schneider and 

G.J. Milburn, |quant-ph/01 12042 , H. Mack and M. Freyberger, 

juant-ph/0206177| ; M.B. Plenio and S.F. Huelga, Phys. Rev. 

Lett. 88, 197901 (2002); S. Furuichi, M. Ohya and H. Suyari, 

Rep. Math. Phys. 44, 81 (1999). 
[2] S. Bose, I. Fuentes-Guridi, PL. Knight and V. Vedral, Phys. 

Rev. Lett. 87, 050401 (2001) and erratum, ibid. 87, 279901 

(2001). 

[3] J. Eisert a nd M.B. Plenio, Phys. Rev. Lett, (in print), q uant- 
ph/01 11016. 

[4] J. Gemmer and G. Mahler, Eur. Phys. J. D 17, 385 (2001). 
[5] E.T. Jaynes and FW. Cummings, Proc. IEEE 51, 89 (1963). 
[6] B.W. Shore and PL. Knight, J. Mod. Opt. 40, 1195 (1993). 
[7] W. Vogel, D.-G. Welsch, and S. Wallentowitz, Quantum Optics, 

An Introduction (Wiley- VCH, Berlin, 2001). 
[8] M. Horodecki, P. Horodecki, and R. Horodecki, Phys. Rev. Lett. 

80, 5239 (1998). 

[9] K. Zyczkowski, P. Horodecki, A. Sanpera, and M. Lewenstein, 
Phys. Rev. A 58, 883 (1998); J. Eisert and M.B. Plenio, J. Mod. 
Opt. 46, 145 (1999). 

[10] G. Vidal and R.F Werner, Phys. Rev. A 65, 032314 (2002). 

[11] J. Eisert, PhD thesis (Potsdam, February 2001). 

[12] R. Clifton, H. Halvorson, and A. Kent, Phys. Rev. A 61, 042101 
(2000). 



[13] P. Horodecki, J.I. Cirac, and M. Lewenstein, juant-ph/0103076 ; 
J. Eisert, C. Simon, and M.B. Plenio, J. Phys. A 39, 3911 
(2002). 



5 



